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Abstract 

We substantiate our statement that the deuteron remains bound in the chiral limit. We critically 
discuss recent claims that effective field theory cannot give a definite answer to this question. 
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1. The question of the quark mass dependence of the nuclear forces has attracted considerable interest 
recently. The tool to perform such investigations is chiral effective field theory (EFT), which allows 
for a unified description of pion-nucleon and nucleon-nucleon dynamics. In Q we had shown that 
the deuteron remains bound in the chiral limit, making use of an EFT approach based on a modified 
Weinberg power counting which successfully describes two-, three- and forur-nucleon systems. In 
earlier 0] and parallel work Q a different claim was made, namely that within the EFT one presently 
cannot make a definite statement whether or not the deuteron is bound in the chiral limit because 
some low-energy constants (LECs) are not known. In these works, a coordinate space power counting 
was used which is supposedly equivalent to the one employed by us. In a more recent paper, Beane 
and Savage || tried to substantiate their claim by performing an error analysis for the various input 
parameters similar to what was done already in JjJ. Here, we wish to critically examine various claims 
made in that paper, and as a consequence, we are able to substantiate our original statement. 

2. Since we do not wish to repeat any of the detailed arguments underlying our calculation presented 
in [P, we rather concentrate here on some issues raised in || which deserve a closer look. 

1) In Ref. ||, the naive dimensional analysis for the quark (pion) mass dependent operators DjM% 
(where / specifies the channel under consideration) is repeated in terms of a parameter rj. It is 
claimed that we use rj < 1/19 in the 3 Si channel. We are puzzled by such a statement. Even 
though the four-nucleon LECs are to some extent regulator and convention dependent, such 
effects should largely cancel in ratios of the LECs projected onto partial waves. So what we have 
really used for this relative strength is 

_ Di So Ml 1 _ Ds Sl Ml 1 

which is more than a factor of two different to what is claimed in [||]. The values of the LECs 
C S ° and C Sl corresponding to the exponential regulator defined in |J and the cut-off A = 560 
used here and in [|J read: 

Co So = -0.143 x 10 4 GeV" 2 , C 3 Q Sl = -0.146 x 10 4 GeV" 2 . (2) 

The appropriate (not partial wave projected) LECs Cs and Ct, which enters the underlying 
chiral Lagrangian Q, take the values: 

C s = -0.979 F~ 2 , C T = -0.003 F~ 2 , (3) 

where we have used the same notation (and normalization) as in Ref. §0 While the con- 



stant Cs is perfectly natural with ac s = —0.979, the LEC Ct appears to be very small as a 
consequence of the approximate Wigner symmetry, see || for more details. The authors of || 
should provide better information on the size of their LECs in the appropriate partial waves. It 
is now important to address the question whether the range for rjj used in |l|] is too narrow? 
As we have explained there, we have performed detailed and serious studies of the size of LECs 
appearing at NLO and NNLO in the chiral EFT and none of the LECs which can be determined 
from data exceeds ±3 when properly normalized in units of Fl and A 2 ,. Therefore, the choice for 
rji made in [0, which corresponds to —3 < a\ < 3, is well rooted in the phenomenology of four- 
nucleon LECs. Also, if one were to extend this range to say ±10, one would leave the realm of 



* 4 Some of the values given here are slightly different from the corresponding ones in || due to the different value of 
g-nN used, see Jl). 
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a converging EFT because NLO correction would become comparable to the LO contributions. 
Thus, even as one claims to know nothing about the quark mass dependent LECs, one should 
still bound them based on such a convergence argument. Finally, as we have discussed in Q, 
approximate Wigner symmetry requires small values of the corresponding LECs. While such an 
argument is merely aesthetic, it still should not be discarded completely. 

2) The central result of Ref. || is Fig. 4. This figure as it is shown is more than confusing. First, all 
other figures are scatter plots, so why are suddenly continuous lines drawn? Simply connecting 
the most outside lying points by a linep 5 ! makes not much sense. In particular, the upper line 
explodes with decreasing pion masses. A simple extrapolation of this upper line to smaller 
values of M n would probably lead to an unnaturally large deuteron binding energy of the order 
of few hundreds MeV and thus indicates the breakdown of the EFT. This could indicate that the 
authors have considered a region of parameter space that does not exclude deep bound states 
outside the range of EFT. The authors of Q do not clearly discuss the origin of such a dramatic 
increase of the deuteron binding energy for smaller values of M % . In Fig. |l|, we show our results 
for the deuteron binding energy taking the range of die as given in the caption of Fig. 4 of M and 
increasing the value of r]3g l given in eq. (|]) and used in Q by a factor of two. The resulting value 
r)3g = 1/4.3 is even somewhat larger than the one suggested in ||. The deuteron is bound in 
the chiral limit and the binding energy only varies between 4 and 20 MeV.pl In our calculation, 



the expansion of the binding energy is well behaved for all values of M n . The scattering length 
in the 3 Si partial wave calculated with the uncertainties in the LECs as suggested || is depicted 
in Fig. H 

3) There is an apparent inconsistency between Figs. 2,3 and Fig. 4 of ||. If for pion masses below 
60 MeV the deuteron can be bound or unbound, the scattering length in the 3 Si channel should 
diverge as the deuteron passes from being bound to being unbound provided the M„--dependence 
of the deuteron binding energy is continous. No such divergence is seen in these figures. We 
also remark that one rather should plot bands (as done in jlj]) because the variations of these 
parameters cannot simply be treated as statistical fluctuations. Furthermore, as discussed in 
0, there are sizeable correlations between certain LECs in that analysis. It is therefore highly 
questionable to use the MINUIT errors as given in that paper for a statistical error analysis as 
done in ||. 

4) It is not clear from Ref. || that the calculation is done consistently with the power counting 
advocated in It that scheme one should treat the difference between the OPE exchange 
and its chiral limit representation perturbatively. This is a very slowly converging expansion for 
small energies as demonstrated explicitly in . As shown in Eq. (3) of || , the full OPE is used 
in H and no uncertainty due to this procedure is given. 

5) Keeping the short-range terms DjM% together with the chiral limit representation of the two- 
pion exchange (TPE) in the potential as is done in H, pi seems to be inconsistent. The reason is 
that the renormalized expression for the leading chiral two-pion exchange contribution includes 
apart from the non-polynomial pieces also terms of the kind [jl|] 

ai M% + frM* ln^, (4) 

#5 As was explained to us by one of the authors. 

* 6 We do not consider a variation in dig as relevant here. The value we use refers to the well accepted small value of 
the pion-nucleon coupling constant as deduced by the Nijmegen and VPI/GWU groups. The larger uncertainty due to 
the fitting procedure in M should therefore be taken cum grano salis. We will come back to that point later on. 
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Figure 1: Deuteron binding energy as a function of the pion mass. The shaded 
areas correspond to the allowed values. The light shaded band refers to the variation 
of D 3Sl using r)s Sl = 1/4.3 and <f 16 = -1.23 GeV~ 2 , d 18 = -0.97 _GeV" 2 . The dark 
shaded band gives the uncertainty if, in addition to variation of D3 g x , the LEC die 
is varied in the range from die = —0.17 GeV -2 to die = —2.61 GeV -2 . The heavy 
dot shows the binding energy for the physical value of the pion mass. 



where oti and Pi are the known constants and A refers to the renormalization scale. ^_ This 
indicates that the M^-dependence of the leading TPE is as important as the M^-dependence 
of the short-range terms in eq. (||). It is, therefore, not clear, why the authors of ||, || decided 
to neglect the explicit M^-dependence of the TPE as well as the second term in eq. (Q) and to 
keep only the first term in that equation. 

6) The authors of || claim to be able to reproduce our results using the same input parameters. As 
we just showed with respect to the discussion of their Fig. 4, we are not able to reproduce theirs. 
In particular, we obtain for the deuteron binding energy in the chiral limit: Bffi = 9.6 ^4 
MeV, where the uncertainties for D3$ 1 (the first error) and die (the second error) are taken 
from |H or are even slightly larger, i.e.: rj3g 1 = 1/4.3; — 2.61 GeV -2 < die < — 0.17GeV -2 . 
As already pointed out, we do not consider a variation in dis as relevant here, and, therefore, 
have not plotted the corresponding bands in Fig. |]. For the sake of completeness, we however 
calculated the resulting additional uncertainty in the chiral limit value of the deuteron binding 
energy. Note that the LEC dig does not contribute to the OPE in the chiral limit (where the 
Goldberger-Treiman relation is exact) and thus changing Jig only affects -B§ L indirectly, due 
to corresponding small changes in the LECs related to contact interactions, as explained in 
Q. Variation of dig in the range —1.54 GeV -2 < dis < —0.51 GeV~ 2 in addition to variation of 



* 7 Similar terms also arise from renormalization of the short-range interactions by pion loops. 



4 



100 200 300 400 500 600 700 



Figure 2: 3 S% -scattering length as a function of the pion mass. The shaded areas 
correspond to the allowed values. Bands as in Fig. [|. The heavy dot shows the 
scattering length for the physical value of the pion mass. The triangles refer to 
lattice QCD results from ||. 



D-js-l and di$ as suggested in [|| leads to Bp L = 9.6^ 3 ' 2 ^ 4 1q 4 MeV. As expected, the resulting 
additional uncertainty is quite small. 

7) Some of parameter choices discussed in H appear questionable. While all fits of die m {§ give 
negative values for this LEC, the authors of || claim that a positive value of d\e = +1 GeV -2 is 
consistent with the data. It is their duty to present an analysis of the many data on ttN — ► ttttN 
using this particular value for d\Q. It should also be stressed that the various fits of ttN — ► ttttN 
in [H are for exactly one choice of irN phases and thus one value for dig (the phase shifts of H ) . 
Therefore, the independent variation of d\e and d\%, the latter referring to other ttN phase shift 
analyses, as done in Ref. || does not really make sense. 

3. In summary, we have shown that the results presented in [jl]] are robust under parameter variations, 
in particular, the deuteron is bound in the chiral limit. We have also shown that some of the results 
presented in || cannot be understood simply and thus need clarification. We stress again that the 
discussions presented here do not influence any of the statements made in our previous work all 
results and conclusions of that paper remain unchanged. 

Acknowledgements We thank Silas Beane for some clarifying discussions. 
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